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1. Introduction 

In this paper we treat the split Stirling refrigerator driven by a 
linear compressor. From the many dissipative mechanisms in the 
system only the friction of the compressor piston, the friction of 
the regenerator, and the viscous losses due to the gas flow through 
the regenerator matrix will be taken into account. The temperature 
at the cold end is taken as an input parameter. First the general 
equations are derived which are subsequently treated in the har¬ 
monic approximation. In principle this limits the application of 
the formalism to the steady state and to low amplitudes. On the 
other hand it has great advantages since it leads to simple analyt¬ 
ical formulae which leads to a clear picture of the basic operation 
of the machine. They allows a quick overview of the essential as¬ 
pects in the design, and demonstrates fundamental parameter 
combinations. In addition the relations have elegant solutions 
which correspond to the operation as a cooler, as expected, but also 
as a heat engine. 

In a split Stirling refrigerator the regenerator/displacer is driven 
by the dynamic pressure in the system. There is a finite pressure 
drop over the regenerator. At the same time a finite flow through 
the regenerator is needed to produce the cooling. The combination 
of the pressure drop and the volume flow is a fundamental source 
of dissipation. In the absence of a bouncing volume the regenerator 
is driven only by the pressure drop over it. In that case the split 
Stirling refrigerator needs a fundamentally dissipative process to 
generate the cooling. In this respect it resembles a pulse-tube 
refrigerator which needs an orifice. This type of cooler is called a 
Stirling refrigerator but the cycle differs from the fundamental Stir¬ 
ling cycle which is a reversible cycle between two isotherms and 
two isochores and has the Carnot coefficient of performance 
(COP). In this paper we will derive a relation for the optimum 
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COP of the split Stirling refrigerator without bouncing volume 
which is fundamentally lower than the Carnot COP. Fortunately 
the degradation of the performance due to the dissipation in the 
regenerator can be made rather small in practice. 

The theoretical analysis of split Stirling cryocoolers has been the 
topic of various papers. The analysis starts with a paper by de Jonge 
[1]. It gives a phenomenological description in which the system 
properties such as pressure variations are related to the positions 
of the piston and the regenerator. The actual values of the phenom¬ 
enological parameters are derived from a numerical computer pro¬ 
gram and from experiment. In this sense the phenomenological 
parameters are used as fitting parameters. Ackermann [2] analyzed 
free-piston resonant cryorefrigerators in the harmonic model 
emphasizing the use of the phasor formalism. Since then many pa¬ 
pers have been published related to this topic. A number of them is 
given in Refs. [3-9]. For designing practical refrigerators our work 
results in some general rules of thumb. For more detailed design 
and optimization more elaborate numerical calculations are 
needed such as in [10] or with Sage and Regen [11]. 

Coolers and heat engines have much in common. As mentioned 
above we will show that, for certain combinations of frequency and 
temperatures, the split Stirling refrigerator actually operates as a 
free-piston Stirling heat engine. A full discussion of the system as 
a heat engine can be found in the book of Walker [12]. 

2. System description 

2.1. Components 

A split Stirling cryocooler consists of a compressor, usually of 
the dual-piston flexure-bearing type, connected to a cold finger 
via a split pipe. Fig. 1 represents the cooler in linear form. The exit 
of the compressor and the entrance of the cold finger both have a 
heat exchanger to ensure that the gas, shuttling between these two 
components through the split pipe (with negligible volume), 
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Nomenclature 



A 

area (m 2 ) 

*0 

regenerator position amplitude (m) 

D 

molar flow conductance (mol/sPa) 

X 

general dynamic system parameter 

F 

force (N) 

z 

specific flow impedance (1/m 2 ) 

f 

friction factor (N s/m) 



K 

axial length parameter (m) 

Greek symbols 

H 

enthalpy flow (W) 

oc 

area ratio 

k 

spring constant (N/m) 

y 

specific heat ratio 

L 

length (m) 

n 

viscosity (sPa) 

m 

mass (kg) 

V 

frequency (Hz) 

N 

help parameter (Eq. (54)) 

S 

coefficient of performance 

n 

molar flow rate (mol/s) 

777 

help parameter (rad 2 /s 2 ) 

P 

pressure (Pa) 

(P 

phase angle (rad) 

P 

power (W) 

CO 

angular frequency (rad/s) 

Q 

quality factor 



Q 

heating or cooling power (W) 

Subscripts 

R 

ideal gas constant (J/Kmol) 

0 

time average 

Si 

entropy production rate (W/K) 

1.2, b, 

c, d, e, f spaces in the system 

t 

time (s) 

c 

Carnot 

T 

temperature (I<) 

cc 

relative to Carnot 

V 

velocity (m/s) 

g 

gas 

¥ 

volume (m 3 ) 

L, H 

low, high temperature 

V 

volume flow rate (m 3 /s) 

r 

regenerator, resonance 

w 

compliance factor (m 3 /Pa) 

Q 

see Eq. (65) 

X 

position (m) 

X, y 

in phase, out off phase 



Fig. 1 . Schematic diagram of the split Stirling cryocooler indicating the sign con¬ 
vention and the labeling of the various components and spaces. 


always is at room temperature. In the cold finger the regenerator 
also acts as the displacer. (These functions can also be separated.) 
We will simply call it the regenerator. The regenerator is fixed to a 
guiding rod and a mechanical spring. The guiding rod can stick into 
a backing volume (displacer gas-spring volume). The mechanical 
spring can be mounted directly behind the regenerator or inside 
the backing volume but that makes no difference for the analysis. 
The motion of the regenerator is driven by the pressure difference 
over it (between spaces d and e) and by the pressure difference 
over the guiding rod (between spaces d and f). The compressor is 
represented as a single mass-spring system driven by an external 
force F applying an average power P. How the power P is generated 
and how the overall coefficient of performance (so including, e.g. 
the electrical losses) can be optimized is an interesting topic. How¬ 
ever, in this paper, we concentrate on the cooler part and consider 
the power as an input parameter. 

Many of the parameters describing the system are defined in 
Fig. 1. In mechanical terms the split Stirling cryocooler is a system 
of two coupled oscillators set in forced oscillation by the external 
force F. One oscillator is the piston with its spring and gas-spring 
action from space b. The other oscillator is the regenerator with 
its spring and the gas-spring action from the spaces d, f, and e. 
The cold head absorbs heat at a rate Ql at a temperature F L . Room 


temperature will be denoted by F H . Most of the time we will as¬ 
sume that F l < F h but by the end of the paper we will see that 
there are special solutions of the relations which correspond to 
heat-engine operation, so we will also indicate what may happen 
if the cold end is heated to temperatures above room temperature 
(F l > F h ). 

The pressures in the system all vary around an average value 
p 0 =p according to p = p 0 + dp. We assume that the pressure vari¬ 
ations are small compared to the average pressure |<5p| <p 0 . We 
also assume that the variations of the various volumes due to the 
motion of the piston and the regenerator are much smaller than 
the corresponding average volumes. These two assumptions allow 
linearization of the relations. We also assume that the pressure in 
each of the volumes is homogeneous and that there is no pressure 
drop over the split pipe so that 

Pd=Pc- (1) 


We will meet several times the situation where the working 
fluid flows into a certain control volume V with volume flow V A 
and flows out with volume flow V 2 while the pressure p inside 
the volume changes with time t. If the gas in the control volume 
is an ideal gas V\ can be related to V 2 by 


Vi 


V dp 
y'p d t 


* 


+ V 2 - 



This important relation is valid even if the temperature distribution 
in the control volume is not homogeneous. A volume flow in Eq. (2) 
can be replaced by a moving wall with surface area A and velocity v. 
In that case we have to replace V by vA. If the process in the volume 
is adiabatic the factor / = y = C p /C v which is the ratio of the heat 
capacities at constant pressure and at constant volume. If the pro¬ 
cess in the control volume is isothermal / = 1. If the process is in 
between isothermal and adiabatic some empirical value between 
1 and 7 can be used. In order to simplify notations we introduce 
the compliance factor w, = V i0 /y'p 0 (with i = b, c, d, e) and the molar 
volumes V H = RT H /p 0 and V L = RT L /p 0 , with R the molar ideal-gas 
constant. 
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2.2. Compressor 

For the backing volume b in the compressor holds 
V b = ^bo + A c x c (3) 

and for the volume c at the compression side 

V c = V c0 - A c x c (4) 

with x c the displacement from the equilibrium position of the pis¬ 
ton and A c the area of the piston. Using the dot notation for the time 
derivative the rate of change of the volumes b and c are related to 
the piston velocity v c = x c by 

\/ b = -\/ c =A c x c . (5) 

It is assumed that there is no significant flow of gas to the back 
space. In that case with the w r notation introduced above 

£14 = A c x c = -w b <Sp b ( 6 ) 

and differentiating 

V b = -w b p b . (7) 


Here A r is the area of the regenerator which is equal to the area A e 
the moving surface of space e which in turn is equal to the sum of 
the areas of space d A d and of the guiding rod Af 

A y = A e = Ad + Af (15) 

and x r the displacement from the equilibrium position of the 
regenerator. 

Eq. (2), applied to volume d, gives 


v d = V c - WdPd - V H n r . 

(16) 

With Eqs. (10) and (12) 


* 

A d x r = Vc - W d p c + V H D r p r - 

(17) 

For volume e 


V e = Vifl Y - W e p e - 

(18) 

From Eqs. (1), (10), (11), and (13) 

A r X r = ViDrPr + W e Pr + W e p c . 

(19) 


2.4. Powers 


This relation is a special case of Eq. (2). Eq. (2), applied to the com¬ 
pression space c, gives 

. * 

V c = -WcPc - V c ( 8 ) 

* 

with V c the volume flow of the gas leaving the compressor. With Eq. 

(5) 

* 

A c x c = w c p c + Vc. (9) 

2.3. Cold finger 

In practice the void-volume fraction (porosity) in the regenera¬ 
tor may be between 30% and 70%. The hydrodynamic and thermal 
properties of regenerators with a nonzero void volume can be dealt 
with by integrating expression (88)-(92) as derived in [10]. Also 
numerical programs such as Sage or Regen can be used [11]. Under 
somewhat idealized conditions even analytical expressions, using 
Bessel functions, can be obtained [14]. However, these expressions 
are rather complicated and far from transparent. Therefore, in this 
treatment, we assume that the void volume in the regenerator is 
zero. This will limit the applicability of our approach. However, 
the influence of the void volume of the regenerator can be included 
to some extend by adding part of the void volume of the regener¬ 
ator to the volumes V d0 and V e0 . In the case that no gas is stored 
inside the regenerator the molar flow n r of the gas flowing through 
the regenerator is the same everywhere in the regenerator (what 
flows in flows out). In the linear (low velocity) approximation n r 
is given by a Darcy type of relation 

n r = -D r p r (10) 

with 

Pr = Pe-Pd (11) 

and D r the molar conductivity of the regenerator. 

The volumes of spaces d, e, and f are given by 

V d = V d0 + A d x r ( 12 ) 

and 

V e = V e0 — A r x r (13) 

( 14 ) 


In this section we derive a relation for the power P in terms of 
the dynamic parameters of the system. We start with the equation 
of motion of the piston, with F the external force, which reads 

m c x c = F + A c (<5p b - Sp c ) - kcXc -f c x c . (20) 

Here k c is the spring constant of the piston suspension (e.g. flexure 
bearings) and/ c the friction coefficient. With Eq. (6) we can write 

F = m c x c + A c 3p c + k' c x c +/ c x c (21) 

with 

4 2 

/<=—+ /<c (22) 

w b 

the effective spring constant. The first term is due to the elasticity 
(gas spring) of the gas in space b. The average power, performed 
by the external force F on the piston, is given by 

P = W C . (23) 

With Eq. (21) 

P = x c x c m c + Sp c v c A c + x c x c k' c +/ c v2. (24) 

In this paper we will use several times the relation 


XX = 0. (25) 

This holds since in the steady state 


dX_ ldX 
dt — 2 dt 


(26) 


Eq. (25) confirms that, on average, no net power is needed for a 
periodic reversible energy conversion. Using Eq. (25) for X = x c 
and X = x c respectively shows that the first and third terms in Eq. 
(24) are zero. So the conversions of the kinetic energy and the 
spring energy require no net power input. With this result and 
Eqs. (5) and (24) can be written as 


P=-Sp c Vc+f c vl (27) 

The first term in Eq. (27) can be expressed in the dynamic 
parameters of the cold finger using Eq. ( 8 ) 

* _ 

P= ( 5p c (w c p c + \/ c )+/ c V 2. 

Using again Eq. (25), now for X = p c , this reduces to 


and 

Vf = Vfo + AfX r . 


( 28 ) 
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P = P w +/ c v2 (29) 

where 

Pw = Sp c V c (30) 

is the so-called pV-work. With Eq. (16) 

Pw = V d Sp c + V H h r 5p c . (31) 

With Eqs. (1), (10), and (11) 

Pw = Vd<5p d + V H h r p e + D r V H pi. (32) 

With Eqs. (10), and (25) (for X = p e ) we get 

Pw = C7^ + (AC7p e + D t V H tf. (33) 

1 L 


The last two terms in Eq. (33) can be recognized as due to the 
dissipation by the flow in the regenerator and due to the friction 
of the piston. The first two terms can be cast into recognizable form 
if we use the expression for the acceleration of the regenerator 

m r x r = p 0 A f + p d A d - p e A r - k' r x r -f r k r . (34) 

Similar to the derivation of the gas spring due to space b the gas in 
space f leads to a gas-spring force on the regenerator, so we 
introduced 

K = + k r . (35) 

1 w f 

Multiplying relation (34) with v r , and taking the time average and 
using v/ = 0 and Eq. (25) for X = x r and X = x r = v r we end up with 

p^r x A d = p~ e V r A r +/rV2 (36) 

or 

p d V d = -pjr e +f r ^. ( 37 ) 

Substitution in Eq. (33) gives 

Pw=g- lVe^Pe + Dr^HP? +/ r V?- (38) 

The cooling power Ql we derive from the first law applied to the 
space e 

0 = Ql + H g e - V e 3p e . (39) 

The last term in the right-hand side of Eq. (39) is the pV-work. The 
second term is the net enthalpy flow of the working fluid flowing 
through the regenerator. If it can be arranged that the fluid enters 
and leaves space e with a constant temperature (T L ) we have 

H ge = 0 (40) 

so the cooling power can be expressed as 

Q l = VJf e . (41) 

Substitution of Eq. (41) in (38) gives 

P w= (^-i)q l + D r V H p2 +/ r v?. (42) 

Eq. (42) with (29) is a special form of the general relation! 13] 

P= ^-i1qi + T h S(. (43) 

The first term is the well-known Carnot contribution. The second 
term is the dissipated energy with S, is the entropy production 
due to irreversible processes in the system. In our case 

(44) 


In fact Eq. (42) could have been written down immediately based on 
the first and second laws of thermodynamics and without the 
rather lengthy derivation given here. The agreement between Eqs. 
(42) and (43) is an important consistency check. 

2.5. Summary 

We will now summarize the results obtained so far. The average 
power P is given by Eq. (23) or Eq. (42) and the cooling power by 
(41). The force F is given by Eq. (21). The molar flow through the 
regenerator is given by Eq. (10). Writing w cd = w c +w d . we get 
form Eqs. (9) and (17) 

A d x r = A c x c - w cd p c + V H D r p r . (45) 

Reminding Eq. (19) 

A r X r = V L D r p T + WePc + W e p r . (46) 

The equation for the motion of the regenerator is given by Eq. (34). 
Using Eqs. (1) and (15) 

A r p r = -m r x r - A f Sp c - k' r x r -f T x r . (47) 

So we have the three Eqs. (45)-(47) for the position of the compres¬ 
sor piston x c , the pressure drop over the regenerator p r , and the 
pressure in the compression space p c . In the next section these three 
relations will be solved within the powerful framework of the har¬ 
monic approximation. 

3. Harmonic approximation 

In this section we will first derive the general relations form the 
relations given above. Next we consider a situation which we will 
call the basic situation which demonstrates the properties of the 
split Stirling cryocooler in its purest form assuming that the area 
of the guiding rod is zero, A f = 0 and the friction coefficients of 
the piston and of the regenerator are zero (f c = 0 and f r = 0). 

3.1. General relations 

We will now derive the in-phase and out-off phase (lower index 
x and y, respectively) components of the various dynamic parame¬ 
ters describing the system. All variables will be decomposed in 
their harmonic components, e.g. the pressure drop p r over the 
regenerator is written as 

p r = p rx cos cot + p ry sin cot (48) 

with co = 2nv the angular frequency and v the frequency. We will 
define the phases of all dynamic parameters with respect to the po¬ 
sition of the regenerator by substituting 

x r = x 0 cos cot (49) 

with x 0 the amplitude of the position of the regenerator. So the in- 
phase and out-off-phase components of the regenerator position 
are given by x rx = x 0 and x ry = 0, respectively. The velocity compo¬ 
nents of the regenerator are v rx = 0 and Vry = -cox 0 . 

Equating the in-phase harmonic components in Eq. (47) results 
in 

A r p rx = -Tu x Xo - A f p cx , (50) 

with the notation m x co? = k' r and w r = m r (co' r 2 - co 2 ). Equating the 
out-off-phase harmonic components in Eq. (47) results in 

ArP r y = -A f p cy +f x cox 0 . (51) 

The components of the compressor pressure can be obtained from 
Eq. (46) resulting in 

-A r coX o = ViDrPjy - w e cop cx - 


T H Si = D r V H pi +/ r v2 +/ cV 2. 


w e cop rx 


(52) 
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and 

0 = V L D r p rx + w e cop cy + w e cop ry . (53) 

Eqs. (50)-(53) form a set of four linear equations forp rx ,p ry ,p cx , and 
p cy . With the notation 

N = (V L D r A f ) 2 + (w e coA d ) 2 , (54) 

and using Eq. (15), the components of the pressure drop over the 
regenerator are given by 

Np^ = -w e co 2 [(A f A r + w e ro r )A d + A f V L D r / r ]x 0 (55) 

and 

Np ry = co[(w e co) 2 f r A d - V L D r Af(A f A r + w e ro r )]x 0 . (56) 


icc = j~- (69) 

4. Simplified cases 

In this Section we will disregard the friction in the piston so that 
P = Pw- Furthermore we consider two important simplified cases: 
the first one is the case where there is no guiding rod. The driving 
force for the motion of the regenerator is the pressure drop p r over 
it. The second special case is where the flow resistance of the 
regenerator is zero. The motion of the regenerator is due to the 
pressure difference over the guiding rod which has a nonzero cross 
section in this case. 


The amplitude of the pressure drop over the regenerator satisfies 



{AfA Y + nj r w e ) + {cof Y w e ) 22 

(V l DA) 2 + (W e oAi ) 2 0 


The lower index A is used to denote the amplitudes. 

The components of p c are given by 

NPcx = [w e m 2 {A d {A 2 + w e ra r ) + A v V L D r f v ) - V 2 L D r m r A { }x 0 

and 

Npcy = m[ViD r A r (AfA r + w e ro r ) + V 2 D r A f /r - (w e G)) 2 / r A d ]x 0 . 


(57) 


(58) 

(59) 


The relations for the position of the piston can be obtained form 
Eq. (45). This results in 

A c coX cx — A^coXq + w cd ft>p cx + Vy{D r p ry (60) 


and 


A c coX cy — w cd <x»p cy Vf\D r p rx . 


(61) 


The components of the force can be obtained from Eq. (21) 

P x, y — [m c a c + A c p c + k c x c + /cV c ] xy . (62) 


The cooling power can be obtained with Eq. (41). Since v rx 
cooling power is 

Ql = -^r(Pcy +PryKy. 

Substitution of the equations given above 



1 A w n Ad(ft7 r w e + A v Af) +Aff v ViD r 2 2 

—A Y ViU r -~^- CO Xn 

2 (V,DA) 2 + (WeOjAi) 2 


0 the 
(63) 


(64) 


Eq. (64) shows that the cooling power is zero at a transition angular 
frequency coq given by 


4.1. No guiding rod 


In this case Af = 0. Hence A d = A r and k' r = k r . With Eqs. (54) and 
(57) we get for the amplitude of the pressure drop over the 
regenerator 



(70) 


Note that the pressure drop over the regenerator does not depend 
on the molar flow conductance D r of the regenerator. With Eq. 
(64) we get for the cooling power 



1 

2 


ViD r m r v2 


(71) 


The cooling power is zero for w r = 0 so if the applied frequency is 
equal to the mechanical resonance frequency of the regenerator. 
With Eqs. (67), (70) and (71) we get 


P 




W e / r ft> 2 \ 

V L D r m r J 



(72) 


Since for Ql we need w Y > 0. We see from Eq. (72) that the COP is 
reduced due to the flow resistance even if the flow resistance of 
the regenerator is minimized (D r —> oo). 


4.2. Basic case 


Now we consider a system in which the area of the guiding rod 
is zero, Af = 0 and that the friction coefficients of the piston and of 
the regenerator both are zero (J c = 0 and f Y = 0). In this situation 
the operation of the cooler is shown in its purest form. Therefore 
we call this the Basic situation. Now the relation (70) simplifies to 



a | 
co H 

1 m r w e 



(65) 


The power can be expressed in harmonic components of the force 
and the piston velocity as 


PrA=J* 0 . (73) 

Eq. (71) remains unchanged. The relation for the power (72) reduces 
to 


P = ^(FxVcx+FyVcy) (66) 

using Eqs. (62), (60), and (61). An other way to determine P is to use 
Eqs. (29) and (42) 

P = - l) Q l + 1 V H D r p 2 A + I/ r v? A + I/ c v 2 a . (67) 

As usual the COP is defined as 

e = Y- ( 68 ) 

The Carnot COP is equal to £ c = T L /(T H - T L ). In order to quantize 
the performance of a cooler it is more appropriate to relate the 
COP to the Carnot COP by defining 


P 




W e 7i7 r 



(74) 


With Eqs. (73) and (68) the COP is given by 

£ = =TV- (75) 

With td y = m r (co( 2 - co 2 ) and w e = V e0 /yp 0 the effective temperature 
T' h is given by 



r H + 


TuV 


H v e0 


yp<A m . 


K‘ 


a 



(76) 


We see that in the basic case both the power P and the cooling 
power Ql are proportional to the flow conductance D r of the 
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regenerator. As a result the COP is independent of the flow conduc¬ 
tance. Keeping in mind that, due to Eq. (71), for a cooler co < co' r the 
effective temperature T' H > T H so £ < f c . The reduction of the effi¬ 
ciency is due to the viscous losses in the regenerator. The losses 
can be minimized by keeping the regenerator as light as possible 
(small m r ), have a small volume of the expansions space V e o, And 
a large regenerator area A r . Finally it is advantageous to have co as 
close as possible to co' Y . However, co cannot be equal to co' Y since that 
would lead to an infinite oscillation amplitude x 0 . 

4.3. Large flow conductance 

Now we turn to the case where the flow conductance of the 
regenerator is very large (D r —> oo) but the friction factor of the 
regenerator is nonzero and the motion of the regenerator is driven 
by the guiding rod (A f ^ 0). With Eqs. (54) and (57) we get for the 
amplitude of the pressure drop over the regenerator 

J {A f A r + ra r W e ) 2 + (to/ r W e ) 2 

PrAoo = ViD^Af ° >X °' ( 77 ' ) 

For the cooling power we get from Eq. (64) in the limit D r —> oo 

Qioc = ^rfrco 2 x 2 0 . (78) 

In the expression of the power, Eq. (67) using (77), the term due to 
the flow in the regenerator disappears. What remains is 

P 0O = (^-1 + <x f W (79) 

with af defined by A f = a f A r . 

5. Results 

The relations derived so far contain a wealth of phenomena. 
Here we will highlight some of them. As an example we treat the 
system mainly with the parameters given in Table 1 which corre¬ 
spond with a small split Stirling refrigerator. The spring constants 
of the piston and the regenerator are calculated from the resonance 
frequencies v with k = m(2nv) 2 and the friction factors from the 
quality factors Q by / = 2nvm/Q. 


Table 1 

Input values of the system parameters, characterizing the standard situation 


System parameter 

Symbol 

Value 

Average backing volume 

^b0 

200 cm 3 

Piston diameter 

D c 

15 mm 

Regenerator diameter 

D r 

9 mm 

Length of regenerator 

L r 

4 cm 

Regenerator displacement amplitude 

x 0 

1 mm 

Length of space d 

tdO 

10 mm 

Length of expansion space e 

t e 0 

10 mm 

Applied power 

P 

60 W 

Piston/spring resonating frequency 

V C 

50 Hz 

Regenerator/spring resonating frequency 

v r 

70 Hz 

Applied cooling power 

Ql 

5 W 

Low temperature 

T L 

60 K 

Room temperature 

Th 

300 K 

Average pressure 

Po 

3 MPa 

Piston mass 

m c 

100 g 

Regenerator mass 

m r 

6.5 g 

Guiding rod area ratio 

oc r 

0.25 

Specific heat ratio 

y 

1.67 

Specific flow resistivity 

Z r 

0.36 pm -2 

Viscosity 

11 

20 pPa s 

Quality factor compressor 

Qc 

2 

Quality factor regenerator 

Qr 

1 


5.1. Basic system 

We start with a description of the results for the Basic situation 
in which there is no guiding rod (A f = 0) and the friction factors of 
the piston and the regenerator are zero (f c =f T = 0). 

Fig. 2 is a plot of the power and cooling power as functions of 
the frequency according to Eqs. (74) and (71). We see that for v be¬ 
low the resonance frequency v r of the regenerator (here chosen as 
70 Hz) the system is a cooler (both Q L > 0 and P > 0). Fig. 3 is a plot 
of £ cc the COP related to the Carnot COP, derived from Eqs. (69) and 
(75). We see that £ cc < 1 since there is a fundamental source of dis¬ 
sipation in the regenerator. However, for the parameters chosen in 
this example, the degradation in the performance is only a few 
percent. 

The phase cp c of the piston position is given by tan cp c = -x cy /x cx 
(please note the minus sign). As the phase of the regenerator posi¬ 
tion is zero by definition (Eq. (49)) the phase cp c is also the phase 
difference between these two moving components. Fig. 4 shows 
that cp c < 0 for v < v r and changes sign at the resonance frequency 
of the regenerator. A negative cp c means that the compression takes 



Fig. 2. Power P and applied heating power Q L according to Eqs. (74) and (71) res¬ 
pectively for x 0 = 1 mm and P L = 60 I<. 



Fig. 3. Plot of £ cc , he COP relative to Carnot, as function of the frequency. Due to the 
dissipation in the regenerator it is always less than 1. 
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Fig. 4. Phase of the piston displacement as a function of the frequency for the basic 
situation. 

place when the fluid is in the warm regions in the system and 
expansion when it is in the expansion space e, as it should for prop¬ 
er cooler operation. 

Coming back to Fig. 2 we see that, at intermediate frequencies 
(70-321 Hz), both Ql < 0 and P > 0. This means that heat is ex¬ 
tracted at the cold end and work is recovered at room temperature. 
In other words the system works as a heat engine between 300 K 
and the low temperature T L = 60 K. For stationary operation the 
60 I< point has to be cooled somehow. For v > 321 HzQl < 0 and 
P > 0 the dissipation in the regenerator is so high that heat has 
to be applied at the cold end and power has to applied to the 
piston. 

In Fig. 5 we consider the case where T L > I H . In contrast to Fig. 2 
we see that the power P and the applied heat Ql differ in sign over 
the whole frequency range. At low frequencies the system works as 
heat engine; in the high-frequency range as a heat pump. 

In this subsection we assumed that the friction of the piston and 
the regenerator was exactly zero. It is interesting to see what hap¬ 
pens if the friction is nonzero, but still rather small. As an example 
we choose the quality factors Q c = Q r = 20. Fig. 6 is a plot of the 
resulting power and cooling power as functions of the frequency 



Fig. 5. Power P and applied heating power Q. L as functions of the frequency 
according to Eqs. (74) and (71) respectively for x 0 = 1 mm and P L = 600 K. 



J Q _I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_I_L_I_I_ 

0 100 200 300 400 

is [Hz] 


Fig. 6. Power P and applied heating power Q. L as functions of the frequency 
according to Eqs. (42) and (64) respectively for x 0 = 1 mm,T L = 60 K, and 
Q c = Q r = 20. The frequency ranges where the system works as a cooler and as a 
heat engine are indicated. 

according to Eqs. (67) and (64). We see that for v < 190 Hz the sys¬ 
tem is a cooler (Q L > 0 and P > 0). In the frequency region of 
208 Hz to 315 Hz Ql < 0 and P < 0, so heat is extracted at the cold 
end and work is recovered at the hot end. In other words: the sys¬ 
tem works as a heat engine between 300 K and 60 K. The fact that 
the low temperature is below room temperature is impractical, but 
the principle remains the same. Heat-engine operation may show 
up in a cooler during warm up, where the heating power has been 
switched off. It has the tendency to keep on running. It should be 
noted that Fig. 6 holds for an amplitude of the regenerator of 
1 mm. If the system is let to run freely as an engine it will set to 
a certain amplitude and frequency which depends strongly on 
the load. These matters have to be analyzed thoroughly. However, 
this is outside the scope of the paper. 

5.2. General system 

Now we turn our attention to the behavior of the system with 
parameters as given in Table 1 so with quality factors of Q c = 2 



Fig. 7. Cooling power and applied power as functions of the frequency for 
x 0 = 1 mm, T l = 60 K, and Q c = 2 and Q. r = 1. 
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and Q. r = 1. Fig. 7 is a plot of the power and cooling power as func¬ 
tions of the frequency according to Eqs. (42) and (64). At frequen¬ 
cies above 193 Hz the system stops to work as a cooler. Fig. 8 is a 
plot of £ cc the COP related to Carnot. It has an optimum near 
140 Hz. Naturally it is zero at the frequency where the cooling 
power is zero. 

Fig. 9 is a plot of the ratio of the piston displacement amplitude 
x cA to the amplitude F A of the force. This shows the characteristic 
features of a resonating system with a low quality factor. 

Fig. 10 is a plot of £ cc as a function of the guiding rod area ratio 
af. This graph shows that the guiding rod not only keeps the regen¬ 
erator in position, but it can also be used to improve efficiency. 

Fig. 11 is a plot of £ C c as a function of the specific resistivity of 
the regenerator material. It shows that, in general, there is an opti¬ 
mum value for z r . This illustrates that some flow resistance is 
needed to set the regenerator in motion. If also the nonideal heat 
exchange between the working fluid and the matrix would be ta¬ 
ken into account an other reason for an optimum value of z r would 
be involved since heat exchange and flow resistance are related. 

In Fig. 7 P > 0 at all frequencies, so the system never operates as 
a heat engine due to the friction of the piston and the regenerator. 
The possibility of producing work is higher if the ’’low” tempera¬ 
ture is increased to values above room temperature since the first 



Fig. 9. Ratio x cA /F cA as function of the frequency. 



Fig. 10. Plot of £ cc as function of the relative cross section of the guiding rod. 



z r Uxnr 2 ] 

Fig. 11. Plot of the COP relative to Carnot as a function of the specific resistivity of 
the regenerator. 



Fig. 12. Demonstration of the heat-engine characteristics if T L = 600 K, in the sta¬ 
ndard case for x 0 = 1 mm. In the low-frequency regime heat is applied at I L and 
power P is recovered. 
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term in Eq. (42) becomes negative. If T L > T H the power and the ap¬ 
plied-heat plot is as in Fig. 12. Even for the high-dissipation values 
of Q c = 2 and Q r = 1 the system works as a heat engine in the low- 
frequency range. In the high-frequency range the system is a heat 
pump. 

6. Conclusion 

In this paper we have derived the basic relations for the split 
Stirling cryocooler in the harmonic approximation assuming that 
the void volume in the regenerator is zero but taking the flow 
resistance and the friction of the piston and regenerator into ac¬ 
count. The relations show a wealth of phenomena. Their relative 
simplicity allows a description of the system performance in com¬ 
prehensive physical terms. Special attention is paid to the basic 
system in which the guiding rod is absent and the friction of the 
piston and the regenerator is assumed to be equal to zero. In this 
case the relations become even simpler and e.g. allow an analysis 
in term of a cooler and a heat engine. 
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